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I. INTRODUCTION 

In general relativity, due to the equivalence principle, 
it is challenging to define energy density associated with 
gravitational field. As a result, people turn to define 
the total mass of an isolated system. For asymptotic¬ 
ally flat spacetime, the Arnowitt-Deser-Misner mass at 
spatial infinity can well describe the energy of a grav¬ 
itational system. For stationary spacetime, a timelike 
Killing vector exists, based on which the Komar mass 
was defined. For nonstationary fields where there is an 
asymptotic time translation symmetry, one can define the 
Bondi energy at null infinity. Since we usually work in 
a nonisolated system, for convenience, several quasilocal 
mass definitions have been constructed mm- If we con¬ 
sider cases with more symmetry, it is obvious that most 
of these definitions become equivalent to each other and 
the Misner-Sharp mass is normally a well-posed one [3|. 
However, here we explore interiors of black holes in which 
the spacetime is dynamical and the Misner-Sharp mass 
function shows unusual features, testing the validity of 
this mass definition at such regions. 

Under external perturbations, the inner horizon of a 
Reissner-Nordstrom black hole can contract, and near 
the inner horizon the Misner-Sharp mass grows rapidly. 
This phenomenon is called mass inflation [4] [5] . These 
arguments were extended to the rotating black hole case 
in Ref. [5]. The main cause of mass inflation is that, 
in the vicinity of the inner horizon, due to strong grav¬ 
ity and repulsion, the source field changes dramatically, 
accumulates a large amount of energy, and then mod¬ 
ifies the geometry significantly. To a large extent, the 
mass function is a geometrical quantity. As a result, the 
mass function diverges. In order to get more informa¬ 
tion, some numerical simulations in more realistic mod- 
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els have been performed in Refs. tana. An apparent 
horizon, a null, weak mass-inflation singularity along the 
Cauchy horizon, and a final, spacelike, central singular¬ 
ity were obtained. Spherical collapse of a charged scalar 
field with regular initial data was investigated rigorously 
in Ref. [15] . 

Besides the importance of studying dynamics near the 
inner horizons of charged and rotating black holes, it can 
be instructive to explore an even simpler case: dynam¬ 
ics near the central singularity of a Schwarzschild black 
hole. In Refs. m on, the asymptotic dynamics near 
the spacelike singularity in spherical scalar collapse was 
studied both analytically and numerically, while mass in¬ 
flation was ignored. In Ref. [18], black hole formation 
was investigated analytically and a lower bound for the 
divergence rate of the Kretschmann scalar was obtained. 
In Ref. [19], numerical simulation of spherical scalar col¬ 
lapse showed that near the central singularity the mass 
function blows up and verified the lower bound for the 
divergence rate of the Kretschmann scalar reported in 
Ref. [18] , while analytic work on mass inflation was ab¬ 
sent. 

The purpose of this paper is to emphasize and extend 
our previous results on spherical scalar collapse reported 
in Ref. [20;. We find that, due to the strong backreaction 
from the scalar field on the geometry, the mass function 
diverges in the vicinity of the central singularity and the 
Kretschmann scalar grows faster than in the Schwarz¬ 
schild geometry. We argue that in collapse the Misner- 
Sharp mass is a locally conserved quantity, not providing 
information on the black hole mass that is measured at 
asymptotically flat regions. 

This paper is organized as follows. In Sec. [TTJ the 
framework on spherical collapse is developed. We explore 
mass inflation near the central singularity in Sec. |III| The 
behavior of the Kretschmann scalar near the central sin¬ 
gularity is studied in Sec.[IV] In Sec.[V] conservation laws 
in dynamical spacetime and the behavior of the mass 
function for other collapses are discussed. 

Throughout the paper, we set G = c = h = 1. 
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II. FRAMEWORK 

Our line of thought to proceed with the calculations 
is not significantly different from the one we followed in 
Ref. 20. We study the effects of a scalar field on the 
geometry in collapse. To do so, we provide the standard 
definition of the energy-momentum tensor for a massless 
scalar field ip 


x = Xb (see puna for more details): 



m^ r = 4,nr z 


V 



T$ ] = VvVv - (1) 

In addition to this, we use the double-null coordinates, 
ds 2 = —4 e~ 2a dudv + r 2 (dd 2 + sin 2 ddcp 2 ) 

( 2 ) 

= e 2a (—dt 2 + dx 2 ) + r 2 {dd 2 + sin 2 ddcp 2 ), 

where a and r depend on coordinates t and a u = (t — 
x)/2, and v = (t + x)/ 2. Thus, we can find the dynamical 
equations just as in Ref. [2D]: 

+ r, xx ) - r 2 + r 2 x = e~ 2a , (3) 


g,r = 47T rip 2 r . 

From Eqs. ©-©> it is enough to use a second-order 
Taylor expansion to determine the values of r, a, and 
ip at t = At. The value of g at t = At is determined 
by the definition of g. The spatial range is x £ [0; 10]. 
Set r = i /> jX = 0 at x = 0. The value of g at x = 0 
and the outer boundary conditions are obtained via ex¬ 
trapolation. The spatial and temporal grid spacings are 
Ax = At = 0.005. The finite difference method and 
leapfrog integration scheme are implemented. The nu¬ 
merical code is second-order convergent. Further details 
are described in Ref. [5D . 


*tt + cr iXX + r '“ r ' xx + 4t r(V> 2 t - 1>%) = 0, (4) 


III. MASS INFLATION NEAR THE CENTRAL 
SINGULARITY 


2 

- ip,tt + i>,xx + -{~r,tip,t + r tX ip x ) = 0. (5) 

r 

Constraints emerge from {uu} and {iw} components of 
the Einstein equations, 


+ 2 r yU a, u + 4tt rip 2 u = 0, 

(6) 

+ 2r + 4-7T rip 2 v = 0. 

(7) 


For numerical stability concerns, we use the constraint 
equation © instead of Eq. Q [21] ■ Defining a new vari¬ 
able g as g = —2cr — ln(—r (J, we can recast Eq. ([6]) as 
the equation of motion for g , 

g yU = 4tt ■ — • ip 2 (8) 

r, u 

We impose r^t = = <?,t = ip,t = 0 and V’( r ) = 

a ■ tanh [(r — ro) 2 ] at t = 0 to be our initial conditions 
with a = 0.1 and r 0 = 5. The Misner-Sharp mass m can 
be read from its original definition in Ref. [3], 

5 %r, = ^K t + ^l-y. (9) 

We require r = m = g = Oat the origin (x = 0, t = 0). 
The parameters r, m, and g on the initial slice (x = 
x, t = 0) are determined by numerically integrating the 
equations below along the line bounded by x = 0 and 


The code for collapse is run and a black hole is formed. 
The apparent horizon of the formed black hole has a ra¬ 
dius of tah ~ 1.8. We plot the apparent horizon and 
singularity curve r(x,t ) = 0 in Figs. [l|a) and 0b). The 
Misner-Sharp mass function along the slices (x = 1 ,t = t) 
and (x = 2.5, t = t) is shown in Fig.jljc), from which one 
can see that the mass function diverges near the central 
singularity. 

As discussed in Ref. [22], in the vicinity of the central 
singularity, the ratios between the spatial and temporal 
derivatives of the quantities r, <r, and ip take similar val¬ 
ues and are determined by the slope of the singularity 
curve. With this feature and our numerical results, the 
Einstein equations and the equation of motion for ip ©>- 
© can be reduced as follows: 


rr : tt 

- ~r% 

(10) 

<x,tt 

+ A-Ki> 2 t , 
r ’ 

(11) 

ip,tt 

- — r t ipf 
r 

(12) 

The asymptotic solutions to Eqs. (10)-(12) 

are [2D] 122 

r ^ 

A(P, 

(13) 

(7 R 

B In C + cr 0 , 

(14) 

ip » 

Cln£, 

(15) 
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Figure 1: Solutions near to the central singularity in neutral scalar collapse, (a) and (b): apparent horizon and singularity 
curve of the formed black hole, (c) Misner-Sharp mass function along the slices (x = l,t = t) and ( x = 2.5, t = t). (d) and 
(e) are results for the slice (x = 1 ,t = t). (d) lnr = aln£ + b, a = 0.49948 ± 0.00004, b = 0.2172 ± 0.0005. a = aln£ + b, 
a = -0.1651 ±0.0001, b = -0.083 ± 0.001. (e) ip = alnC + b, a = -0.18129 ± 0.00003, b = 0.1322 ± 0.0003. lnm = alnr + b, 
a = -1.670 ± 0.002, b = -1.258 ± 0.009. 


where /? ~ 1/2 and B ss /3(1 — (3) — 47rC 2 . £ = t 0 — t, 
and to stands for the coordinate time on the singularity 
curve. In Fig. [2j ( means the segments of AD and BC. 
Numerical and fitting results for r and a, and ip on the 
slice (x = l,t = t) are plotted in Figs. 3 d ) and [Tfe), 
respectively. 

Combining Eqs. @, ( fl3| ), and ( |l4| ), we are able to 


rewrite the mass function as 


m=-[ l +e -(r 2 -r 2 J] 


-(1 - J 2 )A 3 e 2 * 0 


C 


— D 


(16) 


^(1 - J 2 )A 3+2L)2 e 2a ° 


-2 D* 


where J(= r^/r^) is the steepness of the singularity 
curve and D = \Z8nC. As shown in Eq. (16), the mass 


function inflates near the central singularity. Our results 
from numerical integration and fitting for m are plotted 
in Fig.JlJe). Note that in Ref. [19] the divergence of the 
mass parameter m near the central singularity was also 
reported based on numerical simulation, while analytic 
work was absent. 

Furthermore, we examine the differential form of the 
mass function. Using Eq. and Einstein equations, one 


obtains mm 

(V) 

where x a represents the coordinates t and x. The trace of 
the energy-momentum tensor on these indices T^ (=T/+ 
T“) cancels for the case of a massless scalar field ip. Unit 
vectors along the radial, temporal, and spatial directions 
are denoted by f, t, and x, respectively, f is normal to the 
contour lines r = const, as shown in Fig. [2] Combination 
of Eqs. ([I]) and ([T7| yields the gradient of m 


dm * 

Vm = —1 + 

at 

= —4irr 2 ■ 


= —Anr 2 


dm . 
!dx 5 

i 

2 6 


(V’ft +WPx)[r,tt + r, x {-x)\ 

2a W?t + ip%)r,tVl + J 2 {~f) 


(18) 


«^(-f). 

AA ^ ’ 

From Fig. [2] one obtains 

. . At■ Ax 

AX = 


At 


Note that 


a/ (At) 2 + (Ax) 2 VlTl 2 ' 


^L = ^L.^r f Vm 2 

AA At AA ’ V + ’ 
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Figure 2: Spatial vs temporal variations near the central sin¬ 
gularity in spherical collapse. AA and —|a|r are orthogonal 
to the r = ri and r = r 2 constant curves. The value of a is 
such that the magnitude of —|a|r has the length as shown in 
the figure. 


then we have 


Am Am AA , 9 1 0rT , , 9 , 9 . . . 

™,r ~ -z— = -TT'ir- ~ -47tr 2 --e 2 + ip 2 x ). (19) 


A?’ 


AA A r 


Equation (19) gives the effective energy density for a 


scalar field in the neighborhood of the central singularity, 




( 20 ) 


This is consistent with the Misner-Sharp mass defini¬ 
tion pi). 


There are noticeable similarities between the central 
singularity in a Schwarzschild black hole and the inner 
horizon in a Reissner-Nordstrom black hole during mass 
inflation. The Misner-Sharp mass function is expressed 
by the metric components. For static Schwarzschild and 
Reissner-Nordstrom black holes, the Misner-Sharp mass 
function remains constant throughout the whole space- 
time. For the cases of collapse studied in Schwarzschild 
and Reissner-Nordstrom geometries, when the backreac- 
tion of the scalar field is neglected, the geometry does not 
change and the Misner-Sharp mass function stays fixed. 
When the backreaction is taken into account, the geo¬ 
metry is effectively modified by the scalar field and mass 
inflation occurs. 


IV. DIVERGENCE OF THE KRETSCHMANN 
SCALAR NEAR THE CENTRAL SINGULARITY 

In Ref. |T 8 ] , spherical scalar collapse was explored ana¬ 
lytically, and a lower bound for the divergence rate of the 
Kretsclimann scalar K^R^apR^ 01 ^] inside the appar¬ 
ent horizon was obtained, 

K > for r < r AH . (21) 

A numerical verification of this expression was provided 
in Ref. |19| . Here, besides verifying this result, we derive 


approximate analytic expression for the Kretschmann 
scalar near the central singularity. 

In the double-null coordinates § , the Kretschmann 
scalar can be written as 

4e 4cr 

K =- r "" 4CT 


[e a + r (<r tt - a tXX ) 

+ 2 e - 2 cr (r 2 t — r 2 ) + (r 2 t — r 2 x ) 2 


2 r 2 (r 2 


XX 

■ /l 2 / 2 2 1 2 2 2 2 2 2 \ 

+ 4r (r t CT t + r tX cT x - r t a x - r%a >t ) 

+ 4r 2 (r xx r, x cr tX + r tt r t a,t + r, tt r, x cr,x 
+ r, xx r jt cr,t - 2r , xt r tt cr iX - 2r xi r x a it )\- 

We verify the relation © on several slices (x = Xi,t = 
t ) with 0.5<:Ti<4. The results for a sample slice (x = 
2.5, t = t) are plotted in Fig. [3j 


+ r 2 tt - 2 r 2 


t) 


( 22 ) 


In spite of the large expression in (22), this can be 


simplified in the vicinity of the central singularity, where 
the metric is expressed by the Kasner solution [ 22 ], 


ds 2 = —dr 2 + x 2 pi dxf. 


(23) 


The proper time r is obtained via 


r = 


/■ 


>d(: 


3 + 2 D 2 


C 


(24) 


dx i, dx 2 , and dx 3 correspond to dx, dd, and sin6 d(j) 
in the double-null coordinates ([ 2 ]), respectively. p\ = 
(-1 + 2D 2 )/(3 + 2D 2 ), and p 2 = p 3 = 2/(3 + 2D 2 ). q[= 
4D/(3 + 2D 2 )] represents the scalar field contribution, 
if) = q In r. The Kasner exponents satisfy p\ +p 2 +P 3 = 1 
and p\ + P 2 + p\ = 1 — q 2 - In the Kasner metric (23), 
using Eqs. (13) and (24), the Ricci and Kretschmann 


scalars can be expressed as 

r 6 D 2 { _ 3 + D 2 ) ^ 2 


(3 + 2 D 2 ) 2 
3D 2 (—3 + D 2 ) 


Q 


—3—2D 


(25) 


^ 12(16 - 40D 2 + 99D 4 - 26D 6 + 3D 8 ) _ 4 

K = -(3 + 2D2)-- T 

3(16 - 40D 2 + 99D 4 - 26D 6 + 3D 8 ) / r \ - 2 ( 3 + 21 +) 


64 


(3)' 


(26) 


The expressions (25) and (26) are consistent with those 


obtained in Schwarzschild-like coordinates under quasi- 
homogeneity assumption in Ref. [IT]. H(= 16 — 40D 2 + 
99D 4 — 26D 6 + 3D 8 ) is positive for arbitrary values 
of D and has a minimum value of 11.7 at D « 0.47. 
When D = 0, the Schwarzschild metric is recovered [22], 


and Eqs. (25) and (26) respectively lead to R = 0 and 
K = 48m 2 /r°, as expected. With Eq. (16), there is 


32 to 2 


(1-J 2 ) 2 e 4ff0 /r\— 2 ( 3 + 21 +) 


(3)' 


(27) 
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Figure 3: Verification of Eq. (211: 7v>32m 2 /r 6 for r < tah- The sample slice is (x = 2.5, t = t). (a) results for the time range 
t € [0; to] with grid spacings Ax = A t = 0.005, where to is the coordinate time for r = 0. (b) results near r = 0 obtained via 
mesh refinement. In the figure, Ratio=iv/(32m 2 /r 6 ). 


Then K and 32m 2 /r 6 are at the same scale. Equa¬ 
tions (21), (26), and (27) show that the backreaction of 
the scalar held enhances the growth of the Kretschmann 
curvature compared to the Schwarzschild case by a factor 
of r~ iD2 . 


V. DISCUSSIONS 


In collapse, in the vicinities of the central singularity 
of a Schwarzschild black hole and the inner horizons of 
Reissner-Nordstrom and Kerr black holes, the dynamics 
is local. The quantity m is just a parameter which var¬ 
ies at each point, not giving global information on the 
black hole mass. One may interpret this issue via an 
analog to Newtonian gravity. In Newtonian gravity, sup¬ 
pose we want to measure the mass M of a source sphere 
with a pointlike test mass m using a torsion balance. De¬ 
note the distance between the two masses by r. Without 
perturbations, the gravitational force between the two 
masses is simply F = GMm/r 2 . However, if another 
mass 5M passes by the test mass, the gravitational field 
near the test mass can become very strong and local and 
is not able to provide accurate information on the source 
sphere. When the perturbation mass is gone or sticks 
to the source sphere, one can measure the mass of the 
source object accurately again. 

In the case of static stars, the mass function is a glob¬ 
ally defined charge emerging from Killing vectors, which 
work as generating currents. Now we argue that the 
locally defined Misner-Sharp mass cannot be promoted 
to a global quantity during collapse. “Kodama’s mir¬ 
acle” that was revealed for the first time in Ref. [23] 
provides a set of conservation laws required for the study 
of local dynamics in a spherically symmetric spacetime. 
Conserved charges might be found from Kodama’s vec¬ 
tor flow, which is described by K a = —e ab r yb , where 
(a, b ) —> (T, 7') and e ab is the two-dimensional Levi-Civita 


tensor for a curved spacetime. The Misner-Sharp mass 
as defined in Eq. ([£]) is an example of a conserved charge 
from the current J° = G ab K b . The last conserved cur¬ 
rent is equivalent to Eq. for the case of a nonva¬ 
cuum solution in general relativity. We provide a sim¬ 
plified argument to show that it is only in the static 
case when nonlocal currents can be defined. These cur¬ 
rents are the so-called Killing vectors, which are given 
by V„£ m + = 0. It is always possible to repara¬ 

metrize the null generators in a way that allows us to 
write = (1, 0, 0, 0), oriented with the “temporal” 
direction. Note that + V„£ ;J = g^o. Therefore, if 
the reparametrized metric is independent of the “zeroth” 
coordinate, the Killing condition (V„£ M + = 0) is 

fulfilled and £ /J can be called a Killing vector. Hence, it 
is not possible to align the null generators with when 
the spacetime becomes dynamical. 

To study the collapse case, we consider a minor modi¬ 
fication of the usual Schwarzschild coordinates, coherent 
with the discussions in Refs. [23H2S], 


ds 2 



dT 2 + 


dr 2 

^ 2 m 


r 2 dfl 2 



dr 2 + 


dr 2 

2m 


2 dn 2 


where m = m(T,r ), A = A(T,r), and the last redefin¬ 
ition is possible only if A depends on T exclusively, in 
a way that dr/dT = e x . Henceforth, the Kodama vec¬ 
tor K a [= —e ab r^ b = —(1, 0, 0, 0)] is aligned with our 
previous prescription of Nonetheless, m being a 

time-dependent parameter, this is an example in which 
Kodama currents do not satisfy the Killing condition. 
Therefore, all the local charges, Q = f J a dE a , obtained 
by integration on a fixed hypersurface E can be defined 
globally only in the stationary case. Consequently, the 
aforementioned divergence of the mass function must not 
be considered as a nonlocal phenomenon, which would be 
more relevant. However, it is worthwhile to mention that 
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time dependence is not an obstacle in order to find con¬ 
served quantities in the neighborhood of a point where 
the currents are well defined. It is widely known that 
the greatest achievement of the local approach is to hold 
the invariance of these quantities even in dynamical situ¬ 
ations. (See Refs. [231 (Mi for explicit details.) Further¬ 
more, it is not required to have complete knowledge of 
the whole evolution history to provide these quantities. 
Such a peculiarity can also represent a limitation under 
local (time and position) displacements of the mass func¬ 
tion, affecting the goal of a unique mass parameter. In 
fact, as shown in Eq. (16 1, to is a function of r and D. 
As demonstrated in Ref. [22], the parameter D , which 
represents the strength of the scalar field, varies slowly 
along the singularity curve. 

Mass inflation emerges from the backreaction of the 
scalar held on the geometry. A massless scalar held is 
also called “stiff matter” with equation of state p = p, 
where p and p denote the pressure and energy density, 
respectively. In fact, not every type of matter held can 
have such a strong backreaction so as to modify the geo¬ 
metry. In perfect fluid collapse with equation of state 
p = kp where 0 < k < 1, the dynamics near a space¬ 
like singularity remains the same as in vacuum up to the 
second order [13 [27]. Considering spherical symmetry, 
spacetime geometry is given by the Schwarzschild solu¬ 
tion and the Misner-Sharp mass becomes a constant. The 


Lemaitre-Tolman-Bondi dust collapse model is a typical 
example, in which the mass function is well defined and 
hnite all the way up to the spacetime singularity [25] , 
Now we summarize our results. First, besides the 
inner horizons of charged and rotating black holes, it 
was found that mass inflation also takes place in the 
vicinity of the central singularity of a neutral black hole. 
Secondly, it was argued that, for a collapsing system, 
the quasilocal Misner-Sharp mass definition cannot be 
extended to describe the global dynamics inside the 
black hole. Thirdly, we obtained approximate analytic 
expression for the Kretschmann scalar near the central 
singularity and found that the scalar field backreaction 
makes the Kretschmann scalar grow faster than in the 
Schwarzschild geometry. 
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